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^~^ ' Abstract: In this paper, we consider the numerical convergence of G-heat equation which was first 

^vq ' introduced by Peng. The G-heat equation extends the classical heat equation with uncertain volatility. 

\—i • For G-heat equation is nonlinear partial differential equation(PDE), we prove that the Newton iteration is 

■ ■ convergence and the fully implicit discretization is monotone and stable. Then, we have the fully implicit 



< 

in 



discretization convergence to the viscosity solution of a G-heat equation. 
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Z ! 1 Introduction 

jrt ' In recent years, a nonlinear expectation-G-expectation was established by Peng. In the theory of G- 

expectation, the G-normal distribution and G-Brownian motion were introduced and the corresponding 
stochastic calculus of Ito's type were established(see [5], [5], [7], [5], |1]). In Markovian case, the G-expectation 
is associated with fully nonlinear PDEs, which typical applications among economic and financial models 

Qs , with volatility uncertainty (see [3]). 

O^ , Here, we study a kind very interesting nonlinear equation-G-heat equation, which is related to the G- 

normal distribution. G-normal distribution include distribution uncertain, and general the classical normal 

^^ I distribution. In economic and financial market, we usual assume some variables are classical normal distri- 

bution. But in many time it is fall, we need a new distribution to consider more complicated phenomenon. 
G-normal distribution is a new choice. 

Under G-expectation, we relate the G-normal distribution to G-heat equation. Following the work of [3], 

K^ ' [ID]) [H]) [12] 1 we use the technique in [3] to prove that the Newton iteration is convergence and the fully 

implicit discretization is monotone and stable. 

In the last section. Under the same maximum volatility, we compare the value of G-normal distribution 
and classical normal distribution, and find that the values of G-normal distribution is larger more than 
classical normal distribution for the same event. So the G-normal distribution should be more prudent to 
measure risk. The further study see [2]. 
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The paper is organized as follows: in section 2, we give the notations and results on G-expectation. The 
numerical convergence to viscosity solution of fully implicit is established in section 3. Finally, we give an 
example to compare G-heat equation and heat equation in section 4. 

2 Preliminaries 

Firstly, we give the relation between G-hcat equation and G-normal distribution. 

Let J7 be a given set and let H be a vector lattice of real valued functions defined on f2, namely c E H 
for each constant c and \X\ E H ii X E H. H is considered as the space of random variables. 

Definition 2.1 A sublinear expectation M, onT-L is a functional^ : "H — > i? satisfying the following properties: 
for all X,Y E H, we have 

(a) Monotonicity: If X>Y then t[X] > E[Y]; 

(h) Constant preservation: ]E[c] = c; 

(c) Sub-additivity: E[X + Y]<E[X]+ E[Y]; 

(d) Positive homogeneity: E[\X] ~ XE[X] for each A > 0. 

(r2,'H,E) is called a sublinear expectation space. 

Definition 2.2 Let Xi and X2 be two n-dimensional random vectors defined respectively in sublinear ex- 
pectation spaces (Oij'HijIEi) and (5l2j'H2,IE2). They are called identically distributed, denoted by Xi = X2, 
ifEi[ip{Xi)] = E2[(^(X2)], for all (p E Ci,Lip{R'^), where Ci,Lip{R^) is the space of real continuous functions 
defined on R" such that 

Mx) - Viy)\ < C{1 + Ixl'^ + \y\'')\x - y\ for all x,y E i?", 

where k and C depend only on ip. 

Definition 2.3 In a sublinear expectation space (ilj'HjE), a random vector Y = (Yi, • • • ,Yn), Yi E %, is 
said to be independent of another random vector X = {Xi, ■ ■ ■ ,X„i), Xi E H under E[-], denoted by YJ-X, 
if for every test function ip E Ci,Lip{R"^ x i?") we have E[(p{X,Y)] = E[E[(p{x,Y)]x=x]- 

Definition 2.4 (G-normal distribution) A d-dimensional random vector X = (-^^i, • • • ,Xd) in a sublinear 
expectation space (r2,'H,E) is called G-normally distributed if for each a,b>0 we have 

aX + bX i^ v/a2 + b^X, 

where X is an independent copy of X , i.e., X ~ X and X^-X . Here the letter G denotes the function 

G{A):=]^E[{AX,X)]:^d^R, 

where S^ denotes the collection of d x d symmetric matrices. 



Peng [5] showed that X = {Xi, ■ ■ -^Xd) is G-normally distributed if and only if for each Lp e CLLipiR"^), 
u{t, x) := ¥\ip{x + y/tX)], {t, x) e [0, cx)) x R'^, is the solution of the following G-heat equation: 

dtu ~ G{Dl^u) = 0, u(0, x) = ip{x). (2.1) 

The function G(-) : S^ -^ i? is a monotonic, sublinear mapping on S^ and G{A) = ^W,[{AX,X)] < 
||A|]E[|Xp] —: ^\A\a^ implies that there exists a bounded, convex and closed subset F C S^ such that 

G(A) = isuptr[7A], 

where §J denotes the collection of nonnegative elements in S^. 

Let d = 1, we consider the finite difference method to the next G-heat equation: 



which 



and 



dtU-^{a^{Dl^u)+ -a^{Dl^u)-) = 0, x E R, t > 0, 
u{0,x) = ^{x), X £ R, 






(2.2) 



3 Numerical Convergence 

Firstly, we consider a bounded boundary problem of (|2.2p . i.e., 

dtU-\{a^{Dl^u)+-g_\Dl^u)-) = i), 

u{Q,x) = ip{x)^ a;£[a, fe], (3-1) 

u{t, a) = g{t), u{t, b) = h{t) t e [0, T]. 
where (p, g, h are continuous functions. 

3.1 A Finite Difference Discretization 



The equation (j3.ip can be discretized by a standard finite difference method with variable timeweighting to 
give 

<+i - < = ea^K+, 2< + u^,] + (1 - e)a1+'K^l - 2<+i + u^], (3.2) 



where 

CT(r;')^At, 



a- 



2{xiJ^i~Xi){xi—Xi--i) ' 



CT, if rf > 

a(rf):= { , (3.3) 

CT, if rf < 

» ■ (2;i+i-a;i)(a;i-i:i_i) ■ 



In this paper, we consider the fuUy iniphcit schemes with 9 = 0, i.e., 



w»+1-m"= a"+Mu"+' 



^ ..,+, 2<+i+u»+i], (3.4) 

The set of algebraic equation (13. ip is nonhnear for the formula of a{T'^). So we consider the discrete 
equation at each node as 

<^r := < - <+' + «-!,^[^^I^ - 2<+^ + z.:!+^]. 

Following the work of D.M. Pooey [3] (more see. Pang and Qi [TU];Qi and Sun jjj; Sun and Han [H]), 
we must specify the element of the generalized Jacobian that will be used in the Newton iteration. We define 
the derivatives as 

da{TfV ^ I *'. ifr>0 

^^ [^^ ifr<o ' 

For further analysis the Newton iteration, we rewrite the discrete equation p.2p in matrix form. Let 

For convenience, we modify the first and last rows of M as needed to handle the bounded boundary 
conditions. By the discretization schemes in (J3.4I) . the matrix M is a diagonally dominant matric with 
positive diagonals and non-positive off-diagonals. Note that all the elements of the inverse of M are non- 
negative. The discrete equation p.4p can be rewritten as: 

[/ + M"+i]C/"+i =[/", (3.5) 

where / is the identity matrix. Next we prove the convergence of the Newton iteration for full implicit 
schemes. 

3.2 Convergence of the Newtion Iteration Schemes 

For the matrix AI is a diagonally dominant matric, we can analysis the Newton iteration of equation (j3.5p . 
We adopt the Newton timestep as the following scheme: 



(a) Let {U"+^)" = [/"; 

(b) For fc = 0,1,2,- •• Solve 

[/ + M((C/"+i)'=)](C/"+i)'^'+i == [/" (3.6) 

where ([/"+i)'=+i is the {k + l)th iteration, and M{{U''+^)'') means M be dependent on ([/"+!)'=. 

(c) For a given small number e, if 

max |«+i)'=+i - «+i)'=| < e • max(l, |«+')'=+M). 

(d) we end the scheme. 

We give the convergence results about the above Newton iteration as followed: 

Theorem 3.1 The nonlinear iteration liS. 6]J convergence to the unique solution of i3.5\) . for given intial 
iterate {W+^f = C/". 

Proof. For notional convergence, we denote M'^ = M((t/"+^)'^) and U'^ = (t/"+^)'^. So equation p.6p can 
be rewritten as 

[/ + M'=]f7'=+i = [/". (3.7) 

Firstly, we show that the sequence {t/'^}o<fe is monotonically. The k iteration of equation p.7p gives 

[/ + M'=-1][/'^' = [/". (3.8) 

Subtracting equation (|3.7I) from equation p.8p . we have 

[/ + M'']{tj''+^ - U^) = [Af'^'^i - M^]U''. (3.9) 

We consider the right side of l3.91 for each i 



where 



i (xi+x-Xi){xi-Xi-i)^ 

il^ ■.= [u^a,u\,--- ,ui]', 

CT, if ff > 
CT, if ff < 



By the equation ^^, if ff < 0, a{V'yf = ct, then 

At,(a(ff)2-a(f^i)2)^, 



2 
Similarly, if f f > 0, (j{t^,Y = ct, then 

At,(a(ff)2-a(f^i)2)., 



Ff >0; 



Ff >0. 



For the matric / + M*' is a diagonally dominant matric, the inverse of matric / + M^ is non-negative, 
we have 

^fc+i _ ^fc > 0, fc > 1. (3.10) 

Next, we need to prove the sequence {t/'''}o<fc is bounded. Set C,nax = maxiu", Cmin — miniu", 
Umax = maxi u^ , C/min — niiui u^ . By the equation p.8p , we have 

"f-"-''["m-2'iHd]-<, (3-11) 

where 

"j := TTT TT T- (3-12) 

2(a;j+i - Xi){Xi ~ x^^i) 

By the equation p.lip . and Sj^^ > 0, then 



and 



So 



(1 + 2a1~^)u1 < 2at'Um.. + C„ 



and 



9»fc-l 
^i _ fc — 1 i^^'^ ^^ '-^max; 

l + 2a, 

'fc > ^"» /-> . I r* • 
1 + 2d*^^^ 

Set maxi bk^^r = 6i, max^ bkTrr = 62, and < 61, &2 < 1, then we have 

l+2a^ ' l+2a* ' _ ij ^ j 

77 < c' fi ■ > r" • 

'-^inax _ -, , ^--'max5 '-^mm ^ -, , '^min- 
1 - 61 1 - O2 

Now, we prove the uniqueness. Suppose there are two solution to equation p.5p . C/"^^ and U2^^, such 
that 

[/ + Mi]C/f+i = [/", [/ + M2]C/2"+^ = [/". 
Similar the proof of the monotonicity sequence {t/'^}o<fc, we have 

[/ + M2]iU^'+' - U^+') = [Ml - M2]K+\ 

and 

jjn+l _ jjn+1 > Q 

By the equality of U^+^ and U^+\ we have U^+^ - U^+^ = 0. 
Thus, we complete the proof. 



3.3 The Convergence of Fully Implicit 

In the above section, we have proved the convergence of the Newton iteration for the nonhnear equation 
(|3.4p . Next, we would to prove the full implicit schemes convergence to the viscosity solution of (|3.ip . By 
the work of Barles in [T], we know that a stable, consistent, and monotone discretization will convergence to 
the viscosity solution. 

Theorem 3.2 The fully implicit discretization hS. 5\) convergences to the solution of the equation \3.1\) . as 
At, Ax -^ 0. 

We first give some important lemmas for prove Theorem 
Review the discrete equation at each node as 



< := < - ur' + o^tiKti - 2<-^^ + <^r'] (3-13) 



then at each step 



<^r«++i\<+\<+i\<)=0, Vz. (3.14) 

In the case fo nondifferentiable (p^, we use the following definition of monotonicity: 

Definition 3.3 A discretization of the form ^3.14^ is monotone if either 
or 

w n+l n+1 n+1 ^n > n 

Next, we prove the monotonicity of the fully implicit discretization. 
Lemma 3.4 The fully implicit discretization i3.13\) is monotone, indepdent of any choice of At and Ax. 
Proof: For any given e > 0, we just to chek the next two euation: 



By the definition of (^", we have 



and 



< < - <+i + a-^l [u-+l - 2<+i + <+i] 

This completes the proof. 
Proof of Theorem [3721 

By the results of Barles, we just to check that the fully implicit discretization is consistent, stable, 
monotone. Fristly, the formula (I3.14p is a consistent discretization. Then Theorem 13.11 shows that the fully 
implicit discretization is monotone. So we need to prove the discretization is stable. Set 

Ul,, = max(max C/f ), g", /i"), U^,,, = min(min C/f ), g", /i"). 

where 5", /i" is the boundary value of the nth times step. Using the same mathod as in Lemma l3.4[ we have 
the more exact results: 

,n+l 



"-^mm — "i — "^max 



Thus, we complete the proof. 



3.4 The Superlinear Expectation 

For reader convenience, we still use the same notions as in sublinear expectation (G-expectation), and give 
the main results of superlinear expectation. 

Definition 3.5 A superlinear expectation M on H is a functional E : H — > i? satisfying the following 
properties: for all X,Y Cz H, we have 

(a) Monotonicity: If X > Y then E[X] > E[Y]; 

(b) Constant preservation: E[c] = c; 

(c) Sub-additivity: E[X + Y]<E[X]+ E[Y]; 

(d) Positive homogeneity: E[\X] = XE[X] for each A > 0. 

(r2,'H,E) is called a sublinear expectation space. 



The bounded boundary problem is 



dtU-^Ja'{Dl,u)+-aHDl,u)-)^0, 



u{0,x) = ip{x), X G [a, fe], 

u{t, a) = g{t), u{t, b) = h{t) t e [0, T]. 



(3.15) 



where f,g,h are continuous functions. 

The equation p.lSp can be discretized by a standard finite difference method with variable timeweighting 
to give 



,n+l n 



0anut,, - 2< + uf_,] + (1 - e)ar'Kt^' - 2<+^ + «I^], 



where 



^i^f) 



'^{^i+l—^i){^i—Xi — l) ' 

CT, if rf > 

a, if rf < 

« • {Xi^i-Xi){Xi-Xi-i)' 

In this paper, we consider the fully implicit schemes with 6 = 0, i.e. 



(3.16) 



(3.17) 



^r'-<= <+M^.i^-2<+i+u^], 



(3.18) 



The set of algebraic equation (I3.15P is nonlinear for the formula of cr(r"). So we consider the discrete 
equation at each node as 



^r := < - <+i + «^ K+l 2<+i + U-+1] 
For further analysis the Newton iteration, we rewrite the discrete equation p.l6p in matrix form. Let 



[/"+! - K+l,<+\... ,<+!]', (7" = K,<,--- ,< 



[M"C/"], :=-<[<+!- 2< + <_i] 



[I + M"+']U''+' ^U", 



(3.19) 



Theorem 3.6 The fully implicit discretization \3.19]) convergences to the solution of the equation h3.15]) . 
as At, Ax — > 0. 



4 Numerical Example 

In this section, we give an example which is important for financial market. The nonlinear probability 



(i, x) :— ^[Ix^^x<y]' (^' ^) S [0' °°) ^ ^' should be the viscosity solution of the following G-heat equation: 

dtu - G{Dl^u) = 0, u(0, x) = I^<y. (4.1) 

i.e., 



dtU-\{a'^{Dl^u)+ -q^{Dl^u)-)^Q, x e R, t > 0, 



2 

u{0,x) = Ix<y, X e R, 
Next, we consider a boundary problem of (14. 2p . i.e 



(4.2) 



9t«-i((Z?L^^)+-i(i?L«)-)-0, 

u{0,x)^I,<o, xe [-9.99,10], (4.3) 

M(t,-9.99) = 1, w(t, 10) = 0, ie[0,l]. 

For a given probability space (fi, J^, P), the linear probability u(i, x) := -E'[^2,+^x<^], {t, x) G [0, oo) x i?, 
should be the viscosity solution of the following heat equation: 

dtu - Dl^u = 0, w(0, a;) = h<y. (4.4) 

We also consider a boundary problem of (|4.4p . i.e., 

dtu-\{{Dl^u)+-{Dl^u)-)^Q, 

M(0,a;) = /,<o, X €[-9.99,10], (4-5) 

u(t,-9.99) = 1, M(t, 10) = 0, ie[0,l]. 
Comparing the value of u(l, x) and w(l, x), x G [—9.99, 10] : 

The value of u at t=1 
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The red line is the value of function u(l, x), x e [—9.99, 10], and the blue line is the the value of function 
u{l,x), X e [-9.99, 10]. By u{t,x) := W^+^x<y] and u{l,x) = E[I^_^_^^^y], we have 

E[Ix<o] = u(l,0) = 0.6680, P{X < 0) = u(l,0) = 0.5010. 
These results show that the G-expectation have more improtant application in financial market (see [5]) 
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